Size-effects in the Density of States in NS and SNS junctions 
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The quasiparticle local density of states (LDOS) is studied in clean NS and SNS junctions 
with increasing transverse size, from quasi-one-dimensional to three-dimensional. It is shown 
that finite transverse dimensions are related to pronounced effects in the LDOS, such as fast 
oscillations superimposed on the quasiparticle interference oscillations (for NS) and addi- 
tional peaks in the bound state spectrum in the subgap region (for SNS). Also, the validity 
of the Andreev approximation is discussed. It turns out to be an acceptable approximation 
in all situations tested. 



I. INTRODUCTION 



In 1964 Andreev described a new kind of reflection process Ijy which electrons incident on a Normal-Metal Super- 
■ conductor (NS) interface are reflected as holes, and vice versatl. This process, now known as Andreev reflection, led 
shortly after its discovery to both theoretical and experimental work on tunneling transport and the related local 
density of states (LDOS) in small superconducting structures involving at least one NS interface. The diniftasions 
of the samples perpendicular to the current flow were essentially macroscopic and the corresponding theoriea3'0 were 
QQ , three-dimensional (3D). _ 
J — ' Miniaturization of devices led to the development of mesoscopic physicsu. The initial model approaches in that new 
, field were one-dimensional (ID). Attention-was focused on an obviously nonequilibrium property, the conductance, 
' and guided by Landauer's early result for itEnj. Finite transverse dimensions were considered by cowiting the number 
of transverse modes. The extension of Landauer's formula to the NS system was given by Lamberttl. In these studies 
only the total density of states of ID systems aaters, being inversely proportional to the velocity. In calculating the 
Josephson current in SNS. junctions Beenakkem applied a more advanced expression for the total density of states 
, given by Akkermans et aEl. ITLin 

As far as the LDOS is concerned, even recent studiesE2rt3 are 3D as yet. In this paper we calculate the LDOS of NS 
and SNS junctions with finite transverse dimensjons., by this considering effectively ID systems and all jjossibilities 
between ID and 3D. A Green function approachtj'Efl is used, inspired by IshiiEj and Tanaka and Tsukadallj. 
^ \ In case of the NS junction, we first investigate the LDOS in the quasi-lD limit of this junction. "Quasi-ID" means 
O ■ in the limit of transverse system size going to zero. The LDOS is shown to exhibit oscillations as a function of 
\ both energy E and distancft-from the interface. This result reproduces previously observed and analyzed oscillations 
J> . in tunneling experimentsEj'LZI. We then increase the transverse dimensions and fipd the appearance of additional 
oscillations. In progressively refined applications of Scanning Tunneling Microscopyll3 these oscillatory effects in the 
^ I LDOS might well become detectable in the near future. In the case of infinite transverse dimensions the additional 
oscillations disappear again. 

In the SNS junction, we study the LDOS in the normal region for energies below the superconducting gap, and 
find in the quasi-lD junction one bound state. With increasing transverse system size, the number of bound states 
increases. 

A point of discussion in our analysis is the role of the Andreev approximation, which is demonstrated to be a good 
approximation for both the NS and the SNS junctions. Its effect becomes noticeable for large transverse dimensions 
only. 

In section || we give a brief outline of the theory and the model used. The LDOS in NS and SNS junctions is 
discussed in sections ID and IV, followed by some conclusions in section 0. 



II. THEORY 

The Green function method used in this paper is described in Ref.l3, and will be published in a forthcoming 
publicatiorO. We refer the reader to these papers for an extensive description of it, and here only summarize the 
aspects which are of direct importance for the calculation of the local density of states. 
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The Green function describes the various ways of propagation from one point in space r to another one r . Here 
we study clean metallic systems consisting of a few layers, in which scattering only takes place at the interfaces 
between the layers. In the presence of an interface, the total Green function G(r, r ) usually consists of two terms: 
one bulk term, accounting for propagation in the material without any influence from the interface, and a scattering 
contribution froixL-intcraction with the interface. Wp. use the expressions for the homogeneous bulk superconductor 
as given by IshiiEj and follow Koperdraad et aZ.EjO in determining the scattering matrix elements for two simple 
systems: a planar NS junction with only one interface, and a SNS junction containing two interfaces. 

The central quantity of this paper, the local quasiparticle density of states in 3D inhomogeneous superconducting 
strucljires, is calculated from the matrix Green function corresponding to the Bogoliubov equations for quasiparticle 
statesllj. This Green function is the solution of the following matrix equation 

[iiu„To-KT3~Dix)]G{r,r')=S{r~r')To (1) 

where [j„ ~ {2n + Vj^nksT are the Matsubara frequencies, K is the free particle Hamiltonian minus the chemical 
potential /x, in atomic units (?i=2m=l) given by 

and the matrices tq, T3, and D{x) are given by 
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Here A (2;) denotes the superconducting pair potential, which is zero in the normal part of the system. 
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FIG. 1. 3D planar NS junction with finite transverse dimension LyL^ 

We apply the above sketched formalism to systems which are infinitely long in the x-direction, but of finite length 
Ly and in the transverse y- and z-directions, for example to a NS junction with a square transverse cross section 
LyLz as drawn in FIG. |^. Contrary to the translationally invariant situation treated usually0'Ei, in which G(r, r ) 
depends on the differences y~y and z — z only, in the present case the dependence on y, z and the primed coordinates 
is not reduced to differences. This means that variations in the LDOS in the transverse directions survive after taking 
y = y and z = z, the latter substitutions being required in calculating the LDOS. However, such dependence cannot 
be measured and further it depends on the precise preparation of the boundaries. Therefore it is sufficient to take 
the average over the transverse directions. This can be done as follows. First, the boundary conditions are chosen 
such that the Green function vanishes at the boundaries in both transverse directions. Only the functions contribute 
that are proportional to sui(kyy) and sin^k^z), with ky — -j^ and kz = -x^' "■y ^'^'i ""^^ being non-negative integers. 

The expansion of the full Green function in terms of these functions has expansion coefficients G{x, x , ky, ky, kz, k^). 

Subsequently, one puts y ~ y and z — z in this expansion and averages it over the wire's cross section. In this 
way, only terms with ky = ky and k^ = kz survive, the corresponding expansion coefficients of which are denoted 

by G(x,x , ky, kz,iuJn). The variable zw„ is added because the Green function is still a solution of Eq. (^). Finally 
one manipulates Eq. (|l|) according to the expansion and averaging procedure indicated above and one finds that it 
reduces to 

[iujnTo - K^Ts - D{x)] G{x, X , ky, kz,iujn) = 5{x - X )to (2) 



2 



with 



= - 4. and kl^=fi-kl- kl (3) 

Using jproperly normalized functions in the complete sets in the y and z directions, it is found that the quasiparticle 
LDOSO reduces to 

p{x,E) = lim - — — "S^ ImGii(x,x, ky, kz, E + iS) (4) 

Gil is the upper left matrix element of G{x, x , ky, kz, icun) with x — x , and the standard replacement of iw„ by 
E + iS has been appliedtll. Im Gu denotes the imaginary part of Gu and E is the quasiparticle energy measured 
with respect to the Fermi energy /i. Eq. forms the basis of all our present calculations. 

A well-known and often applied approximation in calculations concerning inhomogeneous superconducting struc- 
tures is the so-called Andreev approximation (AA). It was first introduced by Andreevid and can be stated in several 
ways. Perhaps the simplest is to say that in AA normal reflections due to mismatch of wavevectors at the normal 
and superconducting sides of a NS interface are neglected. Mathematically, it means that we make a series expansion 
of the electron and hole wavevectors, and only take into account terms up to first order in E/fi and A/fi. If the 
wavevector appears as a prefactor, it is approximated even further and taken to be the Fermi wavevector. 
In ID systems, where one is often interested in energies E deviating very little from the large Fermi energy, this is 
regarded as a good approximation. In our 3D systems, with an effective chemical potential kp^ = fj. — ky — k^, appli- 
cation of the AA relies on the assumption that A ^ kp^ . It is however not a priori clear whether this assumption 
is valid for all ky and fc^; especially for large transverse dimensions Ly and L^, when p, as given by Eq. (^), is a sum 
over many wavevectors ky, k^, there are terms for which kp is of the same order of magnitude as E and A. 



III. THE NS JUNCTION 

We consider a normal-metal to superconductor junction as in FIG. |l[ From now onwards, in the actual calculations 
the transverse dimensions are taken equal, so Ly = = Lf . In principle, the pair potential A has to be determined 
self-consistentlyt3, but as a first approximation we take it to be zero in the normal region and constant in the super- 
conductor. So the proximity effect is not included. The chemical potentials of the normal metal and superconductor 
are denoted by hn and /is respectively. 
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FIG. 2. LDOS as a function of distance a; in a quasi-lD NS junction with Lt=4; /i]v=/Js=/i=0.5 Ry and A=0.0001 Ry. 
The interface is located at x=0. The inset shows Friedel oscillations in \pN ~ pN,buik\ ■ 10~® as a function of x in the normal 
metal for i?/A=1.01; note that the scales on the y-axis differ by a factor 10^. 



3 



FIG. 1^ shows the quasi-lD LDOS in a NS junction with /iat = /zg = /x, as a function of x for various energies E > A. 
The calculation is exact, i.e. without the AA, and typical values of the chemical potentials and the gap energy in a 
superconductor are used, expressed in atomic units. Oscillations are clearly visible, both in the normal metal LDOS 
{pn, see inset) and in the superconductor LDOS (ps)- The oscillations in p^ are the well-known Friedel oscillations, 
due to interference of incident and reflected electron wave functions, which give rise to a component in pj^ proportional 
to cos(2A:'^a;), where fc"^ — \/ pn + E. The characteristic wavelength of the oscillations is then L^^'' = ~ 4. Since 
there is no potential barrier at the interface, the amplitude of these oscillations is very small, 5 orders of magnitude 
smaller than the ones in the superconductor. These Friedel oscillations would not be fcumd in the AA. 

The oscillations in ps are caused by quasiparticle interference. Mc Millan and RowelO named them the supercon- 
ducting analogue of the Friedel oscillations. Let the electronlike (holelike) quasiparticle wavenumber be denoted by 

(g''); the oscillatory component in ps ~ cos((g'^ — q'^)x) « cos(V E'^ — A?x/ y/ps) then gives rise to a characteristic 
wavelength Lf'^' = J^^., > 2^5^". 



0.05 




FIG. 3. LDOS as a function of E/A in a quasi-lD NS junction at 2: = 3 ■ 10'' in the superconductor. ^j,n=ps=P~0.5 Ry, 
A^O.OOOl Ry and Lt = 4. 

For completeness, we also plot the oscillations in ps as a function oi E/A for a fixed position in the superconductor, 
see FIG. I They were measured in thin films by Rowell and Mc Millantj and TomashllZl. The characteristic energy 

scale is given by = y^^^ + A2. 

Note that despite their "analogous" background, there is also a clear difference between the oscillations in ps and 
the Friedel oscillations. The latter are due to interfering opposite wave vectors of equal magnitude, whereas the former 
are caused by interference of two slightly different parallel wave vectors, and q'^. 

For E < A (not shown in FIG. ) is a decaying function of x, as single quasiparticles cannot propagate into 
the superconductor (evanescent waves). The deca;iM:ate is given by ^^^^'^'^1 VJ^ and the penetration depth is on 
the order of the superconducting coherence lengthc2l, = Ps/{kpA) ^ 0(1000). 
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FIG. 4. LDOS as a function of distance a: in a NS junction for various transverse dimensions Lt of the junction, 
interface is located at x=0. ^=0.5 Ry, A = 0.0001 Ry, _B/A=1.01. 
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Now we increase the transverse dimension Lt- Results in the superconductor for E/A ~ 1.01 are shown in FIG. |^. 
The summation n ™^ LDOS with kp^ = ^ jj, — ky — k1 — ^ ^ — j^{ny + nf) gives rise to fast fluctuations 

superimposed on the above discussed slow oscillations due to quasiparticle interference. As Lt increases, more trans- 
verse modes fit into the transverse dimensions, leading to decreasing size-effects. In the limit — > oo, the summation 
may be replaced by an integration {-^ J^k k ~^ IS^P' / '^^v I which for the LDOS results in (for E > A, and in 
AA) 



pf{x,E)= lira ps{x,E) 

Lt^oo 



= lim 



1 



E 



E-(E- VE^ - A2) cos , ^VE--A 
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where a 

This is the solid line in FIG. ^. It is easy to verify that lima; 
3D bulk superconductor. 



representing the LDOS in a 



In this regime of 3D NS junctions it is interesting to ask what the role of the AA is. As discussed in section ||, 
one would expect this approximation to become worse as the number of transverse modes increases. For the quasi-lD 
system the difference between values of ps with or without AA is typically ~ 10~^%. For systems with small Lt 
{Lt ^ 100) it becomes ~ 0.1% and for large systems at most ^ 1%. This is a factor lO'^ larger than in the ID 
case, although still not visible on the scale of FIG. ^ We are led to conclude that electrons with large transverse 
wavenumbers ky and k^, thus with an angle of incidence deviating considerably from perpendicular to the interface, 
do not contribute much to the LDOS. By using the AA, the value of the LDOS in N reduces and the size-effect 
fluctuations disappear, whereas in the superconductor they are both enhanced. This can be understood by noticing 
that in AA the normal reflections due to mismatch of wavevectors are neglected. On the normal metal side of a NS 
junction with /ijv = /^Si the Fricdcl oscillations, which are caused by normal reflections, are thus suppressed in AA, 
so that pn equals pw.buik- On the other hand, the oscillations in the LDOS in the superconductor, which arc induced 
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by Andreev reflection, are enhanced in AA, due to increased quasiparticle transmission. However, since the amount 
of normal reflection is very small if there is no potential barrier at the interface, the enhancement of ps is also very 
small. 

If there is a potential barrier at the interface, then the dominant normal reflection mechanism is of course not the 
mismatch of wavevectors in N and S, but the presence of the barrier. In that case we would expect that application 
of the AA, even for large Lt, does not lead to significant changes in the value of the LDOS, at most ^ 10^^%. 



IV. THE BOUND STATE SPECTRUM OF A SNS JUNCTION 



E 



N 



FIG. 5. The planar SNS junction: the energy E and pair potential A are measured from the Fermi energy. 



We study a SNS junction as shown in FIG. ^ with the length L of the normal region on the order of the supercon- 
ducting coherence length ^o- The junction has again cross section Lf. The magnitude of the pair potential in both 
superconductors is taken eqpal, but there is a difference in phase 5(j) = (pR — (f>L ■ For a review of this type of weak 
links, we refer to Likharevtll. 

The above Josephson junction is considered with p.^ = p-s = ^- and first in AA. So both interfaces are perfect, and 
there is full Andreev reflection of all quasiparticlcs with i? < A. 
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FIG. 6. Bound state spectrum in the normal part of a quasi-lD SNS junction, calculated within AA. I/t=4, /i=0.5 Ry, 
A=0.0001 Ry, L=2000. 

In FIG. 1^ the LDOS in the normal region is plotted vs. E/A, for various phase differences Scj) and in the quasi-lD 
limit. The continuous spectra for E > A, which are proportional to {fi +_E)^^/^ in the one dimensional case, are 
not shopn,-.Bound states appear at energies satisfying the relation of Kulil£3, which was later also derived by other 
authorsBEa for a ID SNS junction in AA: 

E 

27rn = 2arccos(i;/A) L±5(l> (6) 

In our formalisniEi, Eq. (^) defines the poles of the scattering matrix elements in the matrix Green function; it can 

be understood in the following simple way. 

An electron travelling from x = to x ^ L acquires the phase 

k'^L ~ (j)R ~ arccos(i?/A). 

The first term is the phase accumulated during propagation through the normal metal; the second one is the phase 
shift acquired upon Andreev reflection into a hole and is equal to the phase of the pair potential in the superconductor 
on the right, and the third term stems from evanescently entering of the wavefunction into this superconductorEj. 
Similarly, the back travelling hole acquires the phase 

-k'^L + 4)L~ arccos(£:/A). 

For constructive interference, the total phase acquired on one roundtrip should be an integer multiple of 27r, so 

{k'' ~ k'')L - 2arccos(S/A) + 0l - '/'i? = 27rn 
In AA, so up to terms of first order in -y= in the expansion of fc^ and fc'*, this is the same relation as Eq. (^. 
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FIG. 7. Bound state spectrum in the normal part of a 3D SNS junction, calculated within AA. Lt—S, fi=0.5 Ry, A=0.0001 
Ry, I/=2000 and S(j)—n/4:. The inset shows the same, but then for a large system with Lt— 400. 

Upon increasing Lt one obtains a picture like FIG. 0, which is for Lf — 8, corresponding to two transverse modes. 
In addition to the ID bound state, additional peaks appear due to the finite transverse dimensions. The width of the 
peaks in the figure is determined by the small imaginary part S in the energy E+i(5. The actual bound state energy is 
found for S —^ 0. The number of peaks is equal to the number of different combinations of transverse modes in the y- 
and z-directions. For Lt = 8, the three peaks, from right to left, correspond to the modes with {ny,nz) equal to (0,0) 
(the ID bound state), (1,0) or (0,1) (the bound state which corresponds to one mode in either the y- or z-direction) 
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and (1,1) (the bound state corresponding to one mode in both the y- and z-direction). 

The inset of FIG. shows that for large Lt the dis<:rete peaks due to the finite transverse size of the junction disappear. 
In the hmit Lt oo, one indeed expects a bandp. 

All the above was done in AA. Releasing this approximation yields a bound state spectrum with a slight shift of 
the peaks, as compared to the same calculation in AA. Even for large junctions, this shift is < 10^^ % towards lower 
energies and we thus conclude that the AA is good. 



In conlusion, we have calculated the LDOS in clean mesoscopic superconducting NS and SNS structures with finite 
transverse dimensions. Going from quasi- ID to 3D systems by increasing the transverse dimensions has pronounced 
effects on the LDOS in both types of junctions; in NS junctions additional oscillations are superimposed on the usual 
slow Friedel-like oscillations due to quasiparticle interference. In SNS junctions, wc found additional peaks appearing 
in the bound state spectrum as a function of the transverse system size. 

Besides, we have tested the influence of applying the Andreev approximation, by performing all calculations of the 
LDOS both with and without AA. It turns out, that both in the case of a single NS interface and in the case of a 
SNS junction the AA does not have a large effect on the LDOS, although the AA-induced error grows by a factor of 
10'^ upon going from quasi-lD to 3D systems. It produces a small correction to the value of the LDOS. 

Finally it is worth noting that the Green function method used is not limited to either studies of the LDOS or to 
the mesoscopic junctions considered here; it can also be used to study eg. supercurrents and quasiparticle currents, 
and it can be applied to much larger systems, such as superconducting superlattices. The latter systems were up to 



We would like to thank Mr. P.W. Brouwer, Dr. K. M. Frahm and Mr. J. A. Melsen for critical comments, leading 
to a revised version of the manuscript. Stimulating discussions with Prof. C. van Haesendonck and Dr. A. A. Golubov 
are also gratefully acknowledged. Part of this work was supported by the Stichting voor Fundamenteel Onderzoek 
der Materie (FOM), which is financially supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek 



^ A.F. Andreev, Zh. Eksp.Teor.Fiz. 46, 1823 (1964) [Sov. Phys. JETP 19, 1228 (1964)]; 51, 1510 (1966) [24, 1019 (1967)] 
^ W.L. McMillan, Phys. Rev. 175, 559 (1968). 
3 C. Ishii, Prog. Theor. Phys. 47, 1464 (1972). 

* For a review, see C.W.J. Beenakker and H. van Houten, Solid State Physics 44, 1 (1991). 
^ R. Landauer, IBM J. Res. Dev. 1, 223 (1957). 

® M. Biittiker, Y. Imry, R. Landauer and S. Pinhas, Phys. Rev. B 31, 6207 (1985). 
^C.J. Lambert, J. Phys. Cond. Matter 3, 7331 (1991). 

* C.W.J. Beenakker, Phys. Rev. Lett. 67, 3836 (1991). 

^ E. Akkermans, A. Auerbach, J.E. Avron and B. Shapiro, Phys. Rev. Lett. 66, 76 (1991). 
A. Furusaki and M. Tsukada, Phys. Rev. B 43, 10164 (1991). 

" Y. Tanaka and M. Tsukada, Phys. Rev. B 44, 7578 (1991). 

G. Kieselmann, Phys. Rev. B 35, 6762 (1987); H. Plehn, O.-J. Wacker and R. Kiimmel, Phys. Rev. B 49, 12140 (1994). 

^■^ R.T.W. Koperdraad, The Proximity Effect in Superconducting Metallic Multilayers, PhD thesis Vrije Universiteit, Amster- 
dam, 1995, available on request. 
R.T.W. Koperdraad, A. Lodder, in preparation. 
C. Ishii, Prog. Theor. Phys. 44, 1525 (1970). 

J.M. Rowell and W.L. McMillan, Phys. Rev. Lett. 16, 453 (1966); W.L. McMillan and J.M. Rowell in Superconductivity Vol. 
I, edited by R.D. Parks (Dekker, New York, 1969). 
"W.J. Tomasch, Phys. Rev. Lett. 16, 16 (1966); 15, 672 (1965). 

M.F. Crommie, C. P. Lutz and D. M. Eigler, Nature 363, 524 (1993). 



V. CONCLUSIONS 




ACKNOWLEDGMENTS 



(NWO). 



8 



P.G. dc Gonnos, Superconductivity of Metals and Alloys, Benjamin, New York, 1966. 
M. Kurd and G. Wendin, Phys. Rev. B 49, 15258 (1994). 
K.K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

I.O. Kulik, Zh. Eksp. Teor. Fiz. 57, 1745 (1969) [Sov. Phys. JETP 30, 944 (1970)]. 
H. van Houten and C.W.J. Beenakker, Physica B 175, 187 (1991). 
C.W.J. Beenakker, Phys. Rev. B 46, 12841 (1992). 



9 



